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Abstract. We investigate the elementary excitations of charge and spin degrees for 
the ID interacting two-component Bose and Fermi gases by means of the discrete Bethc 
ansatz equations. Analytic results in the limiting cases of strong and weak interactions 
are derived, where the Bosons are treated in the repulsive and the fermions in the 
strongly attractive regime. We confirm and complement results obtained previously 
from the Bethe ansatz equations in the thermodynamic limit. 
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1. Introduction 

Recent achievements in trapping quantum gases of ultra-cold atoms have opened up 
many exciting possibilities for the experimental investigation of quantum effects in 
low-dimensional many-body systems. In particular, the experimental realization of the 
Tonks-Girardeau gas [1, 2, 3, 4, 5], collective spin wave excitations in the spinor Bose gas 
[6, 7, 8], the quasi ID Fermi gas [9, 10] and the observation of atom matter waves [11] 
have shed light on understanding the quantum nature of ID many-body systems. There 
has been a corresponding revival of interest in the exactly solved models of interacting 
Bose and Fermi gases, such as the Lieb-Liniger Bose gas [12], the Gaudin-Yang Fermi 
gas [13], the ID Hubbard model [14, 15] and the BCS pairing model [16], due to the 
fact that their ground state and low lying excitations are accessible from exact Bethe 
ansatz solutions. 

Quantum gases with multi-spin (pseudo-spin) states exhibit even richer quantum 
effects than their single component counterparts. The update facilities allow one to 
create multi-component Bose gases in low dimensions, e.g., two atomic hyperfine states 
can make up a pseudo-spin doublet [6, 7, 8]. The observation of collective dynamics 
of spin waves in the spinor Bose gas has stimulated theoretical attention on the ID 
integrable spinor Bose gas [17, 18, 19, 20]. In this paper, we present a systematic 
way to derive dispersion relations of low-lying elementary excitations of the spinor 
Bose and Fermi gases via the discrete Bethe ansatz equations in the strong and weak 
coupling limits. This approach has some distinct advantages compared to the usual 
way of deriving the low-lying excitation spectrum: (i) it leads to explicit forms of the 
quasimomenta which are helpful in understanding the collective behaviour of low-lying 
excitations; (ii) it avoids the counting of quantum numbers involved in the integral 
equation formulation of the Bethe equations, and (hi) physical quantities such as the 
charge velocity, spin velocity, effective mass, Fermi momentum, finite-size corrections, 
ground state energy and excitation energies readily follow in the weak- and strong- 
coupling limits. Note that especially in weak coupling limit, those quantites are difficult 
to deal with in the framework of the Bethe ansatz equations in the thermodynamic limit 
(due to singular integral kernels in the linear integral equations). 

To describe elementary excitations in terms of their dispersion relations, consider 
the general case of M. branches (or rapidities) of excitations labeled by fi = 1, . . . ,M. 
For example, Ai — 1 for the spinless Bose gas (charge excitations only), M. = 2 for 
the spinor Bose and spin-| Fermi gas (spin and charge excitations). Generally, the 
relativistic dispersion relation for each one of these excitations is written as 



where A^ is the excitation gap and is the velocity. For vanishing gap A^ = 0, the 
spectrum is linear, such that 



In this case, the underlying effective field theory is conformally invariant. Examples are 




(1) 



(2) 
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the charge channel for spinless bosons with repulsive interaction [21], the charge channel 
for fermions and the spin channel for repulsive fermions. In the other extreme A M ^> 1, 
one ends up with 

v 2 v 2 r> 2 

^HA„ + ^A„ + £, (3, 

which is the classical dispersion of free particles with an effective mass m* . It is known 
that this scenario is realized for spin excitations in the spinor Bose gas with repulsive 
interaction [22, 23]. This case has been analyzed recently [18, 20]. We will consider this 
model with our method in Section 2, thereby confirming quantitatively the following 
results. In the weak coupling limit, the effective mass of the spin wave is almost the 
same as the mass of one boson, while for strong coupling, it tends to the total mass of 
all bosons. On the other hand, charge excitations in this system are still gapless. 

In Section 3, we will turn our attention to the intermediate regime with finite A M , 
which is encountered for spin excitations in the attractive spin-| Fermi gas. In this 
system, the gap is an increasing function of c 2 . The charge excitations again remain 
gapless. For the sake of comparison, we briefly recall the calculation of the velocities for 
the attractive spin-| Fermi gas in the thermodynamic limit at the end of Section 3. We 
note that the method we propose to calculate e(p) in (1) can also be adapted to other 
models such as the ID integrable Hubbard model [14, 15] and the mixed Bose-Fermi 
model [24, 25, 26, 27]. 

2. The spinor Bose gas 

Before turning to the spinor Bose gas, we give an outline of the method we apply. 
2.1. Outline of the method 

In the framework of the Bethe ansatz solution, the energy eigenvalues E and momenta 
p of the two-component Bose and Fermi gases are given in terms of two sets of Bethe 
roots, {kj} . =1 N and {A a } a=1 M , corresponding to charge and spin degrees of freedom. 
Namely, 

N 

E=J2k] (4) 

j'=i 

TV 

p = E k i ( 5 ) 

where {kj} j=1 N and {A a } a=1 M satisfy a set of N + M coupled Bethe equations, 
given below. The ground state momentum po and the ground state energy Eq are given 
uniquely by the location and the number of roots; especially, po = 0. Starting from this 
ground state configuration, low-lying excitations are constructed by changing both the 
location and the number of the roots, thereby obtaining the corresponding momentum 
p and energy E, both parametrized by the root sets. From these, one obtains the 
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dispersion relation E(p) and the excitation energy e(p) = E(p) — E , which is analyzed 
according to Eq. (1). The essential point is that the roots can be obtained explicitly in 
the limits of weak and strong coupling. The thermodynamic limit, which is essential for 
Eq. (1) to be defined at all, is carried out at the very end. 



2.2. Derivation of results 
The Hamiltonian 

i=l 1 l<i<j<N 

describes a quantum gas consisting of N particles with equal masses in one dimension, 
where the particles interact via a pairwise ^-function potential. According to the 
symmetry of the wave function, these may be bosons, fermions or a mixture of both. 
Here, we will consider a two-component Bose gas [18] constrained by periodic boundary 
conditions to a line of length L. The interaction is attractive for c < and repulsive 
for c > 0. The mass is denoted by m and g is the coupling constant which is written 
in terms of the scattering strength c = 2/oid as g — h 2 c/m. An effective ID scattering 
length a 1D can be expressed through the 3D scattering length for the bosons confined 
in a quasi ID geometry. For convenience in calculations, we prefer to set h = 2m = 1, 
and thus to work with the parameter c. However, units may be restored at any time by 
using the dimensionless coupling constant 

7 = mg/(h 2 n) = c/n. (7) 

The energy is given by Eq. (4), where the quasimomenta kj satisfy the Bethe equations 
(BE) [18, 28, 33] 

kj — ki + i c rr K i ~ *u — 2 1C 



t=l J 1 u=l 3 



\ u + he 



/i=l J >~ 2 

1 r — ke — Tjic _ _TT ~ ~ 1 c 
l\ K-h + \ic ~ Xp + ic 

Here j — 1, . . . , TV and a — 1, . . . , M, with M the number of spin-down bosons. 

The ground state is the fully polarized state (total spin S = N/2) with 2S + 1-fold 
degeneracy due to su(2) symmetry [17, 33]. This is achieved by choosing M = in 
Eq. (8). So in the ground state, the BE reduce to those of the Lieb-Liniger Bose gas, 
which have been extensively studied via various methods [29, 30, 31, 32, 34]. From 
previous work [12, 31, 32, 34] on the Lieb-Liniger Bose gas, the ground state energy for 
weak coupling, namely Lc 1, is given by 

E Q w iV(iV - l)c/L. (9) 

In the strong coupling limit, i.e., as the strength of the repulsive interaction tends 
to infinity, the bosons behave like impenetrable fermions and the system is known 
as the Tonks-Girardeau gas. In this regime, the ground state energy E 
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(N 2 — 1) (1 — An/c) follows directly from the asymptotic solutions to the BE [32]. 
Due to conformal invariance, the finite-size corrections to the ground state energy in 
the thermodynamic limit for strong coupling are given by 

E (L, N) - Le^ = + o(l/ L 2 ) (10) 

where v c = vf (1 — An/c) with the Fermi velocity vf = 2nn and the central charge 
C = 1. We notice that in the zero coupling regime c — > the finite-size corrections for 
the ground state energy do not meet the relation (10), because v c = 2y/cn [12] (cf. also 
the discussion after Eq. (13)). This is due to the fact that for c = 0, the Bose gas has 
dispersion h 2 k 2 / (2m), as can be seen from the Hamiltonian (6). In that case the model is 
not conformally invariant (conformal invariance requires a linear dispersion). However, 
for any finite c, one identifies a finite Fermi momentum hp- Then the dispersion relation 
can again be linearized around zbkp, yielding an effective bosonic theory which describes 
the density fluctuations with wave vectors k ~ ±k F [21]. This theory is conformally 
invariant due to the linearized spectrum. 

Spin charge separation is a typical phenomenon in interacting many-body systems. 
We first concentrate on the elementary excitations in the charge sector. These have 
linear dispersion, where the sound velocity is obtained from Eq. (2). In terms of the 
Bethe roots, these excitations are created by particle-hole pairs with respect to the 
ground state configuration of the {kj}- =1 M . Let us first consider the weak-coupling 
limit, such that ki oc y/c/L, i — 1, . . . , N [32]. Then, to order c 2 , the equations 

N 2 N N 

- 1 - 2 E Jjr^ - * E (< E - h ! kt) - h Z ke) (") 

hold, which determine the roots of the BE (8) asymptotically. The lowest possible 
excitations are obtained from the set of equations [32] 

2c ^ 1 



i k N - kg. 



2c N 



kj ~ — V r (13) 

J L ^ i kj -k t 

where j = 1,...,N — 1 and the momentum p parametrizes the excitation. Here, 
p = 2£ir/L with £ integer; however, we will treat p as a variable in the following, which 
will take continuous values in the thermodynamic limit. Note that p is introduced such 
that Eq. (5) holds. The excited energy in the long wavelength limit is easily calculated 
from equations (13), i.e., E(p) m N(N — l)c/L + pks +p 2 , where the Fermi boundary 
ks ~ is obtained from the Hermite polynomial description [35]. Therefore 

E{p) — Eq ps 2y/cnp = v c p. Here the charge velocity is given by v c = 2^/nc, which 

coincides with the result obtained from (9) according to the formula v c = J"^ [|^] 
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[12]. From the Bethe roots (13), the collective effect is that if one particle jumps out of 
the Dirac sea the rest also rearranges. 

Similarly, for the strong coupling limit Lc^> 1, if the largest quasimomentum k^ 
jumps out of the Dirac sea, asymptotic roots of the BE (8) can be obtained from a 
l/(Lc) expansion as 

2iV\ 
Tc~) 



k 



N 



(N-l)ir 2p 



1 



A. j 



2rijTi 2p 
L + Tc 



1 - 



2N 

Tc 



(14) 



where nj = ±1, ±3, 



state energy E 



31? 



. , ±(N-3)/2, -(N-l)/2 if N is odd. It follows that E(p)-E w 
v c p with the charge velocity v c = ^jf- (l — ^) and the ground 
N(N 2 — 1)(1 — ^). Note that as c — > oo the velocity v c becomes 
the Fermi velocity of noninteracting spinless fermions. 

We now focus on the spin excitation with regard to one spin flipping. This 
corresponds to introducing one A-solution into Eqs. (8), which we call Ai. In the long 
wavelength limit a quadratic dispersion of spin wave excitations above the ferromagnetic 
ground state is obtained [22, 23, 20] 

E(p) -E pa p 2 /2m* (15) 

where m* is an effective mass. This is the case indicated in Eq. (3). Now the small- 
coupling expansion of the BE (8) yields 

N 



2c i c c , 

k i ~ ~l 2^ l~. - u. + TT" + Tvi k i 



ke LX 1 L\{ 



(16) 



for j — 1, . . . , N. From Eq. (5), one obtains Ai « ^, which in the long- wavelength limit 
diverges as 1/p, so that Ai 3> k B . We now calculate the energy by inserting Eq. (16) 
into Eq. (4), 

N 
3=1 



E{p) 



L N ^-^ + Lk P 



1 + 



E + p 2 



(17) 



which indicates m/m* pa 1, where m is the boson mass. In the above equation E is 
given by Eq. (9). We notice that higher order corrections to the ground state energy are 
very hard to derive from our expansion ansatz. The next correction to the ground state 



energy of the Lieb-Liniger Bose gas is known to be given by Eq 
Thus from the relation (17) we observe that 

E(p) ^ Eq 



N 2 c 



1 - 



3tt 



1 + 



[31]. 



(18) 



This suggests the effective mass m/m* 



1 - 



3vr • 
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The divergence Ai ~ 1/p continues to hold in the strong coupling regime c ^> 1. In 
this regime, the consistency conditions derived from the BE (8) in the l/(Lc) expansion 
are 

sin ^- 2 £^ + iL( 1 + !). (19) 

Here for convenience we consider N as an odd number. Then 

'l-Ktij 2p c \ f 2N c 

~L~ + Tc + T\ l ) V ~Tc~ + T\\ 

where rij = ±1, ±3, . . . , ±(N - l)/2. From Eq. (5), we find Ai « Nc/Lp > c > kj 
After some algebra we obtain 

N 



B( P ) = E«?«'^^-i)(i-S 

j=l \ / 



, , 1 27r 2 iV / 2iV\ 4 \ / 4iV\ , x 

+ ^h + -3u{ 1+ u) + Li){ 1 -u) (21) 



which suggests the effective mass 

m 1 27r 2 n / 2n 



,1 • (22) 

m* N 3c \ c J y ' 

The leading correction term coincides with the earlier result [20]. The next correction 
term is relevant for finite c. In the strong coupling limit, i.e., c — > oo, all the bosons 
become indistinguishable thus behaving like hardcore bosons with Fermi pressure-like 
kinetic energy. The effective mass takes the maximum value m* = Nm, meaning that 
by moving one boson with down spin, one has to move all the particles with up spins. 
In this way the spin velocity is very small, as if the spin-flip is frozen [20]. 



3. The Fermi gas 

The spin-| Fermi gas with equal number of spin- up and spin-down fermions described 
by the Hamiltonian (6) has been extensively studied in the context of the BCS-BEC 
crossover [36, 37, 38, 39]. Very recently, a comprehensive study of the model with 
arbitrary polarization was undertaken [40]. In order to compare different collective 
dynamics of elementary excitations in the charge and spin degrees with those of the 
spinor Bose gas, we derive the charge and spin velocities for the strongly attractive 
Fermi gas. This is done from an expansion of the discrete BE, similarly to the previous 
section. At the end of this section, we will compare the results with those obtained from 
the BE in the thermodynamic limit. 



3.1. Velocities from the discrete Bethe ansatz equations 

In the strongly attractive regime the bound states behave like bosonic tightly bound 
molecular dimers. The ground state consists of two-body bound states spread out 
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about the origin in pseudomomentum space due to Fermi pressure, i.e., with kj = 
Xj + \\c, k_j = Xj — |ic for j = 1, . . . , N/2, where the parameters {A_,-} . =1 M satisfy 
the BE [13, 30] 

M \ \ _L_ • 

e i2A J L = _TT A J -A, + ic ^ M = N/2 (23) 

for j = 1, . . . , M. The lowest charge excitation comes from the configuration in which 
a bound state with the largest momentum k max = ( M ~ 1 ) 7r (\ _ which is at the edge 
of the Dirac sea in quasimomentum space, jumps out of the Dirac sea. This leads to 
a bound state density fluctuation. The total real momentum changes from zero to p. 
Taking M even, we may obtain the roots of the BE (23) as 



A 



(M-l)n + p + p 



x 3 



2L 2 2Lc 

2rij — 1 p 



M 



2L + 2Lc 



where j = 1, . . . , M - 1 and rtj = — (M - 2)/2, — (M - 4)/2, . . . , (M — 2)/2. It follows 
that 

w-^^X 1 -™)^}-™)* (25) 

It is seen that to this order E(p) — E = v c p with charge velocity 

«. = T (l - -) (26) 
which coincides with the result of Ref. [36] and with the result derived from the ground 
state energy E [36, 37, 39, 40], 

* = - 1 4? N+ m? N i"-v( 1 - ? £)- (27) 

Again we may find from the ground state energy that the finite-size corrections in the 
strongly repulsive limit for the Fermi gas also meet the relation (10) with central charge 
C — 1. The spin-channel does not contribute here, because it is gapped (cf. below). 

We now focus on the spin velocity which involves breaking one bound state over 
the ground state in order to get spin flipping. The ground state is a spin singlet. The 
lowest spin excitation is incurred by a spin triplet excitation with total spin M z = 1. 
In such a configuration the BE become [30] 

M ^ \ x , : 2 , , . i. 



g i2A i L _ _ TT A i ~ A l + 10 TT Aj fc M + 2 ic 
M A,- - X e - ic Xj - hp - \ic 



with j = 1, . . . M — 1 and /i = 1, 2. Withought loss of generality, we assume M is even. 
The BE (28) admit the asymptotic roots 



fci 



n + ir 
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A, 




i^ j + ? (^M Ui _m^v (29) 



In the above equations = 0, ±1, • • • , ±(M — 2)/2 and n± are odd numbers to be 
determined. The number of bound states is now M — 1. This leads to a different parity 
for the BE (28). We can see this directly from the roots (29) where all bound states 
shift by nearly n/2 from their original positions in the ground state configuration. The 
unpaired two fermions are scattered away. Summing up both sides of the third equation 
in Eq. (29) gives the relation Y^=i A? ~ + ^)- This relation implies that 

Yl!f=i \ i s negligible in comparison with ki + k 2 in the strong coupling regime. Further, 
from Eq. (5) we have 

fci + fc 2 ~p(l- 4( ^~ ^ (30) 

if we consider n_ = —n + for minimizing the energy. Using the BE in (28) gives another 
relation 

,-,^-i^). (31 ) 

Thus we find all roots to be given by 

3 \ L Lc Lc ! ) V Lc 

^U^^f.^Mzll). (32) 



L 2/ V Lc 

Here rij = 0, ±1, . . . , ±(M — 2)/2. The quantum numbers n-t for the unpaired fermions 
should be chosen to be as small as possible, i.e., ±1, in order to minimize the energy 
[30]. Actually, n± may take any two odd integer values, which correspond to states 
with energy higher than the ground state energy. After some algebra, we obtain the 
excitation energy 

E(p) « - 2(M - 4 + ^ (l - 2 -^) N ( N-< - 4) 

The choice of n± affects the excitation energy only to the indicated order. From Eq. (33) 
we deduce the dispersion relation e(p) = E(p) — E = A + \ p 2 (l — ^) , where the spin 
excitation gap to leading order in c is A as c 2 /2 and E is given by (27). On the other 
hand, by comparing with Eq. (1), the spin velocity follows as 



(34) 



which is divergent due to the gapped spin triplet excitation [41, 40, 42]. This result has 
been obtained previously in Ref. [36]. 
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3.2. Velocities from the Bethe ansatz equations in the thermodynamic limit 

Spin and charge velocites for the spin-| Fermi gas in the attractive regime have been 
calculated [36, 40], using the BE in the thermodynamic limit. Within this approach, the 
ideas sketched in Section 2.1 are applied to the log form of Eqs. (8), thereby converting 
products to sums. By introducing densities, these can in turn be converted to integrals. 
One ends up with linear integral equations which can be solved numerically in order to 
determine v c and v s . The advantage of this approach is that it is valid for arbitrary 
couplings. In Fig. 1, v c and v s are depicted from a numerical solution of those integral 
equations for attractive fermions. Furthermore, the results (26), (34), are confirmed. 




Figure 1. a) Charge and spin velocity divided by vf = for the attractive regime, 
dependent on 7. The dashed lines are the small 7 approximations, v c . s /vp — lztj/n 2 . 
b) The corrections to v CtS in the strongly attractive limit, namely ~'j(2v c /vp — 1) 
and (it\/2v s /(vfj) + 1)7/2. According to Eqs. (26), (34) and (7), these expressions 
approach 1 in the limit 7 — -> —00, which is confirmed numerically. 



4. Concluding remarks 

In conclusion, we have derived explicitly the dispersion relations for elementary charge 
and spin excitations for the two-component Bose and Fermi gases in the long wavelength 
limit. Analytic results have been found for the charge and spin velocities and for the 
effective mass connected with spin wave excitations in the repulsive spinor Bose gas. 
These quantities reveal significant insight into the nature of the collective hydrodynamics 
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of the Bose and Fermi gases. In particular, physical quantities such as the charge 
velocity, spin velocity, effective mass, Fermi momentum, finite-size corrections, ground 
state and excitation energies, are obtained via an analysis of the discrete Bethe ansatz 
solutions, as opposed to using the continuous integral equation approach. As such we 
believe that it can be easily adapted to a number of related models. 
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